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Abstract. 

We describe the methodology required for estimation of photometric estimates of metallicity 
based on the SDSS gri passbands, which can be used to probe the properties of main-sequence 
stars beyond ~ 10 kpc, complementing studies of nearby stars from more metallicity-sensitive 
color indices that involve the u passband. As a first application of this approach, we determine 
photometric metal abundance estimates for individual main-sequence stars in the Virgo Over- 
density, which covers almost 1000 deg 2 on the sky, based on a calibration of the metallicity 
sensitivity of stellar isochrones in the gri filter passbands using field stars with well-determined 
spectroscopic metal abundances. Despite the low precision of the method for individual stars, 
internal errors of <7[ Fc / H ] ~ 0.1 dex can be achieved for bulk stellar populations. The global 
metal abundance of the Virgo Overdensity determined in this way is ([Fe/H]) = —2.0 ± 0.1 (in- 
ternal) ±0.5 (systematic), from photometric measurements of 0.7 million stars with heliocentric 
distances from ~ 10 kpc to ~ 20 kpc. A preliminary metallicity map, based on results for 2.9 
million stars in the northern SDSS DR-7 footprint, exhibits a shift to lower metallicities as one 
proceeds from the inner- to the outer-halo population, consistent with recent interpretation of 
the kinematics of local samples of stars with spectroscopically available metallicity estimates 
and full space motions. 
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1. Introduction 

The structure, chemistry, and kinematics of the stellar halo of the Galaxy, with its 
predominantly old and metal-poor populations, collectively preserve a detailed record 
of our Galaxy's formation in the early universe. Thanks to large-area surveys such as 
the Sloan Digital Sky Survey (SDSS; York et al. 2000), recent studies have revealed 
that the halo is marked by numerous stellar substructures. The presence of these lumpy 
and complex substructures (both in real space and in phase space) are in qualitative 
agreement with models for the formation of the stellar halo through the hierarchical 
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merging and accretion of low-mass sub-halos (e.g., Bullock & Johnston 2005). Among 
the various substructures discovered to date, the Virgo Overdensity (VOD) is one of 
the most striking. It was discovered as a stellar overdensity of main-sequence stars in 
SDSS; star counts in the region toward Virgo are enhanced by a factor of two above the 
background stellar distribution (Juric et al. 2008). The overdensity seems to be associated 
with clumps of RR Lyrae stars and turn-off stars, but it is less likely to be connected 
with the leading tidal tail of the Sagittarius dwarf galaxy (e.g., Newberg et al. 2007, and 
references therein). At present, metallicity estimation from broadband photometry is the 
only practical means of obtaining metal abundances for a large number of faint objects 
such as those in the VOD. Such methods are based on the relative sensitivity of stellar 
colors to photospheric abundances over a wide wavelength baseline. The clear advantage 
of using a photometric metallicity technique is the efficiency of estimating mctallicities for 
individual main-sequence stars, which are the most plentiful and representative sample 
of stellar populations. 

Ivezic et al. (2008a) constructed photometric metallicity relations in the u — g vs. g — r 
plane using SDSS filter passbands, and studied the abundance structures of the Galaxy 
with an accuracy of ~ 0.2 dex at g < 17. This approach is similar to the traditional UBV 
method (e.g., Carney 1979), which relies on the strong dependence of ?7-band magnitudes 
on metal abundance. However, the u-band photometry in SDSS is limited to u « 22 (99% 
detection limit). This, and the greatly deteriorating errors in SDSS w-band magnitudes 
near the faint limit, restricts photometric metallicity estimates to stars with r < 20.8, an 
insufficient depth to fully explore the VOD. Here we overcome the limitations of the it- 
band photometry in SDSS by exploring less metallicity-sensitive, but better-determined, 
color indices in the gri passbands. A full description of this method and our results can 
be found in An et al. (2009b). Here we only discuss the highlights. 



2. Photometric Metallicity Estimates 

The set of stellar isochrones in An et al. (2009a) were used to determine photometric 
metallicity estimates based on color-color relations for main-sequence stars. A photo- 
metric metallicity ([Fe/H] p hot) was computed for each star by requiring distances from 
main-sequence fitting to be the same from two different CMDs, with g — r and g — i as 
color indices and r as a luminosity index. We searched the entire [Fe/H] grid from —3.0 
to +0.4 to find [Fe/H] with a minimum x 2 value, where we defined 

2 _ (fig-r - jl) 2 (fi g _j - fi) 2 

for each star. Here, [i and are the distance modulus and its error on each CMD, 
respectively. The p, is a weighted average distance modulus from the (g — r, r) and (g — i, r) 
CMDs. In principle, this method can be used for multiple color indices. However, since 
only three passbands are considered here, the problem is reduced to the traditional 
manner of determining metal abundances from a color-color diagram. In some of the 
cases, a star becomes bluer than the main-sequence turnoff as the metallicity increases, 
either due to a large photometric error or a younger age than our assumed values in the 
models. If the minimum % 2 was not found, [Fe/H] p hot was estimated by extrapolating 
the difference in distance modulus toward higher [Fe/H] values based on theoretical 
predictions. 

Although the accuracy of an individual photometric metallicity estimate for a star 
obtained with this technique is rather large (0.5 dex or more), when the method is 
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Figure 1. Left: Number density for over 740,000 stars detected in the SDSS, located at distances 
from ~ 10 kpc to ~ 20 kpc from the Sun in the Lambert projection of Galactic coordinates. 
Each pixel has an area of 12.96 deg 2 , and have a median occupancy of 544 stars per pixel. 
The North Galactic Pole is at the center, and the Galactic Center is at the bottom. Concentric 
circles represent 6 = 0°, 30°, and 60°, respectively. The Virgo Overdensity is the feature seen at 
(/, b) ~ (300°, 70°), while the long strip that crosses the sky from (I, b) ~ (210°, 30°) to the VOD 
is the leading tidal tail of Sgr. The feature at (l,b) ~ (40°, 40°) is the Hercules- Aquila Cloud, 
another large-area overdensity of halo stars discovered in SDSS (Belokurov et al. 2007). Right: 
Median metallicity of the same stars as in the Left panel. 

applied to stellar populations, random errors are beaten down by the large numbers of 
stars involved. In the VOD, for instance (see Figure 1), each pixel in the map is occupied 
by over 500 stars. The median metallicity of stars in the VOD area at 270° ^ I ^ 330° 
and 60° ^ b ^ 70° is [Fe/H] = —2.0 ± 0.1 from the metallicity map in Figure 1, where 
the error is from a pixel-to-pixel dispersion. 

Field halo stars in the mirrored position of the map shown in Figure 1 at 30° ^ I ^ 90° 
and 60° ^ b 70° exhibit ([Fe/H]} = -1.9 ± 0.1. Although half of the stars in the 
direction toward Virgo are likely associated with a progenitor dwarf galaxy or a tidal 
stream, the mean [Fe/H] value essentially remains unchanged; it is the same as that for 
the field halo stars within the precision of the technique. It is perhaps of interest that 
the field star metallicities are as low as they appear to be, as Carollo et al. (2007) have 
argued that the peak of the metallicity of the outer-halo population is [Fe/H]= —2.2, 
and that this component is expected to dominate over the more metal-rich inner-halo 
population (with a peak metallicity at [Fe/H]= —1.6) at Galactocentric distances greater 
than 15-20 kpc. 

3. A Metallicity Map of the Halo to 20 kpc 

The same techniques described above can be used to construct a metallicity map, 
based on the gri passbands, for main-sequence stars in the SDSS footprint extending 
well beyond the limit of 9 kpc reached by the metallicity map presented by Ivezic et al. 
(2008a), based on the ugr passbands. A preliminary result of this exercise is shown in 
Figure 2. We are still in the process of more fully investigating the effect of photometric 
errors at low metallicity, however, the current figure certainly suggests that the typical 
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Figure 2. Metallicity map, obtained from photometric metallicity estimates for 2.9 million 
main-sequence stars with 0.3 < g — r < 0.4 in the northern SDSS DR-7 footprint. The median 
occupancy for pixels shown in the map is 638 stars per pixel. Note that the median metallicity 
for stars located beyond 10 kpc appears much lower than those within 10 kpc. 

stellar metallicity inferred for stars at larger distances is substantially lower than that 
for more nearby stars. 

Carollo et al. (2009) have recently used a new, much larger sample of calibration 
stars from SDSS DR-7 to develop a more quantitative picture of the inner/outer halo 
dichotomy. Their derived velocity ellipsoids for these two components - Inner Halo: (ov R , 
a v<p , a Vz ) = (150 ± 2, 100 ± 2, 85 ± 1) km s" 1 , Outer Halo: (ct Vr , a v<p , a Vz ) = (159 ± 
4, 165 ± 9, 116 ± 3) km s _1 , clearly demonstrate that outer-halo orbits are significantly 
less radially elongated than inner-halo orbits, and possess higher energies, consistent with 
the segregation in spatial dimensions observed in the metallicity map. 

Future imaging surveys, such as Pan-STARRS (Kaiser et al. 2002) and LSST (Ivezic 
et al. 2008b) will use similar photometric bandpasses as those in SDSS, providing even 
deeper (and far more accurate) photometric data than SDSS over a larger fraction of 
the sky. Our photometric metallicity method will be useful to exploit these databases for 
understanding the chemical evolution of progenitor dwarf galaxies that are identified, as 
well as for the bulk populations of field stars. 
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